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SUMMARY 

Contact stresses in meshing spur gear teeth are examined. The analysis is based upon an incremental 
finite element procedure that simultaneously determines the stresses in and the contact region between 
the meshing teeth. The teeth themselves are modeled by two-dimensional plane strain elements. Friction 
effects are included, with the friction forces assumed to obey Coulomb’s law. The analysis also assumes 
that the displacements are small and that the tooth materials are linearly elastic. The analysis procedure 
is validated by comparing its results with those for the classical case of two contacting semicylinders 
obtained from the Hertz method. Agreement is excellent. 


INTRODUCTION 

The life and performance of gear teeth are directly related to the ability of the teeth to withstand 
contact stresses. Contact stresses may produce pitting within the contact area and eventually lead to 
tooth failure. In spite of the importance of contact stresses in gears, comprehensive analyses of these 
stresses have not been extensively reported in the literature. Indeed, most analyses are based upon 
procedures that require simplified assumptions about the geometry of the contacting surfaces. Although 
these assumptions are needed for the classical procedures, their use raises questions about the accuracy 
and applicability of the results. 

In this report we attempt to obtain a more representative and hence a more accurate analysis. The 
method is based upon an incremental finite element procedure (an iterative technique) that simultane- 
ously determines the stresses and the contact area. Our motivation for using the finite element method 
(FEM) is based upon the success of the method in determining fillet stresses due to gear tooth bending 
(refs. 1 and 2). Our motivation for using the incremental procedure is based upon the success of other 
analysts in using it with contact analyses. For example, in 1970 Wilson and Parsons (ref. 3) used a dif- 
ferential displacement method to study frictionless contact problems. This approach was later extended 
by Ohte (ref. 4) to include frictional effects. In 1979 Okamoto and Nakazawa (ref. 5) presented the incre- 
mentation technique. This technique was also presented at about the same time by Skinner and Streiner 
(ref. 6) and by Urzua et al. (ref. 7). The technique was presented in automated form by Torstenfelt in 
1984 (ref. 8). 
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The incremental procedure enjoys a number of advantages over more traditional finite element pro- 
cedures: Specifically, (1) no special elements, such as gap elements or contact elements, are needed; 

(2) the calculated incremental loads follow the actual load history; (3) friction forces , once computed, are 
treated as known tangential loads and thus the standard analysis procedures of the finite element method 
remain intact; and (4) general-purpose finite element codes may readily be used in the analysis. In the 
research described herein we apply the incremental finite element procedure to the special geometry of 
involute spur gear teeth. The balance of this report presents the basic formulation of the method, the 
automation of the method, a validation of the method, some results for spur gear contact stresses, and a 
discussion and some concluding remarks. 
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SYMBOLS 

semilength of contact 

elastic modulus 

accuracy parameter 

nodal force in contact region 

column vector of nodal global forces 

incremental global force vector 

identity submatrix 

global stiffness matrix 

normal unit vector in contact region 

pressure as function of x 

maximum normal contact pressure 

nodal reaction force in contact region 

incremental nodal load vector 

nonsquare transformation matrices 

transformation matrix 

tangential unit vector in contact region 

nodal displacement in contact region 

column vector of nodal global displacements 

incremental global displacement vector 

scale factors defined in equations (17), (18), and (20) 

gap distance 

nodal separation in contact region 
coefficient of friction between contacting bodies 
Poisson’s ratio 
normal stresses 
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r shear stress 

Subscripts and symbols: 
a,b refers to bodies a and b 
i,i' typical pair of contacting nodes 

max maximum 

n,t normal and tangential directions 

represents matrix in local (n,t) coordinate system 


METHOD FORMULATION 


The incremental procedure is an application and natural extension of the basic finite element method: 
Let the governing matrix equation for a finite element analysis be written as 


{F} = [K] {U} 


( 1 ) 


where {F} is the column vector of nodal global forces, [K] is the global stiffness matrix, and {U} is 
the column vector of nodal global displacements. 

Consider a structure discretized into elements and nodes such that equation (1) relates the nodal 
global forces and displacements. Let there be an increment in the nodal forces. Let these incremental 
forces be designated by the column vector {AF}, and let the resulting incremental nodal global dis- 
placements be designated by the vector {AU}. If the global geometry of the structure is essentially 
unchanged by the force increment, the global stiffness matrix [K] is essentially unchanged. Hence, from 
equation (1) the governing equation relating {AF} and {AU} is 


{AF} = [K] {AU} 


( 2 ) 


To apply equation (2), consider two typical bodies A and B in contact over a surface C-C' as 
depicted in figure 1. Let A and B be discretized into elements and nodes in the usual manner of finite 
element analysis, except that in the contact region let each node on the surface of A have a correspond- 
ing opposite, or mating, node on the surface of B (see fig. 1). 

Consider a typical pair of mating, contacting nodes, i and i', as depicted in figure 2. Let t-n be 
local Cartesian axes defining the “average” tangential and normal directions. Let {AR* } and {ARj t } 
(j = a,b) represent incremental normal and tangential loadings on i and i'. Similarly, let {AUj n } 
and {AUj t } (j = a,b) be the resulting incremental displacements. 

The contact status can be categorized as being either (1) open, (2) closed and sticking, or (3) closed 
and sliding. In each case the equilibrium and continuity conditions must be satisfied. That is, 
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(1) For open nodes (no contact) 


{AU b „} - {AU„} < {{„} 

{AR„} - {AR bn } - {0} (3) 

{AR*} - {AR b J - {0} 

(2) For closed and sticking nodes (no relative movement) 

{AU bn } = {AU an } ♦ {6 n } 

{AU at } = {AU bt } 

( 4 ) 

{A Ran) = -{AR bn } 

{AR at } = -{AR b J 


(3) For closed and sliding nodes (tangential movement) 

{AU bn } - {AU„> * {«„} 

{AR bn = -{AR„} (5) 

{AR bl > - -{AR.J - ±^{AR„> 

where {<„> is the normal “gap” vector between contacting nodes and fi is the coefficient of friction 
between the surfaces. (Coulomb friction is assumed and the sign in the final terms of equation (4) is 
chosen so that energy is dissipated.) 

To use the constraint conditions of equations (3) to (5) with the governing equations, it is useful to 
transform the columns and rows, which are associated with the contacting node pairs, from the global 
(x,y) system to the local (n,t) system. Let {AF} and {AU} represent the incremental load and 
displacement vectors in the local system, and let [T] be the (orthogonal) transformation matrix from the 
local to the global system. Then 


{AF} = [T] {AF} 
{AU} = [T] {AU} 


( 6 ) 
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By substituting equations (6) into equation (2), we obtain the relation 


{AF} = [T] t [K][T]{AU} = [K] {AU} ( 7 ) 

where the stiffness matrix [K] is defined by inspection. 

By assembling together the terms associated with the contacting nodes, equation (7) can be expressed 
as 


AF a 


A u a 

AF a„ + AR^ 


A U an 

AF at + AR at 

= [K] 

A U at 

AF bn + AR bn 

A U bn 

AF bt + AR bt 


A U bt 

Xi 

<1 


AU b 


Equation (8) is then the matrix form of the force-displacement relations. To obtain the complete set 
of governing equations, it is necessary to impose the contacting boundary constraints of equations (3) 
to (5) for the three contact conditions: open nodes (no contact), closed and sticking nodes (no relative 
movement), and closed and sliding nodes (tangential movement). 

Open nodes (no contact) . — In this case no constraints are needed. However, equations (3) need to be 
satisfied during each numerical iteration for which any matching node pairs are deemed to be open. 

Closed and sticking nodes (no relative movement) . — In this case equations (4) are constraining rela- 
tions that need to be satisfied. After substitution from equations (4), equations (8) take the form 


AF a 


A u a 

AF an + AR an 


A U M 

AF at + AR at 

= [K] 

A U at 

AF bn - AR an 

A U» - 

AF bt - AR at 


A U at 

AF b 


Au b 


Equation (9) can be reduced by eliminating the unknown contact forces {AR an } and {AR at }. To this 
end, we introduce rectangular transformation matrices [SJ and [S 2 ] defined as 
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( 10 ) 


where I a is an identity matrix with dimensions equal to the number of degrees of freedom of the nodes 
of body a less the number of degrees of freedom of the contact nodes, I b is the corresponding identity 
matrix for body b, and I c is an identity matrix with dimensions equal to the number of degrees of free- 
dom of the contact nodes. For example, if a finite element model has 40 nodes for each of bodies a and 
b and if there are 5 contact node pairs, then I a has dimensions (40 — 5) x 2 = 70, I b has dimensions 
(40 — 5) x 2 = 70, and I c has dimensions 5 X 2 = 10. Observe that these matrices have properties such 
that 


AU a 


AUan 


AU a 


0 

AU at 

AU an ♦ S n 

= [SJ 

AUjui 

+ (^2 ] 


AU at 

0 

AU at 


AU b. 


0 


AU bt 




AF a 

AF a 


AF an ♦ AR m 

AF m+ AF bn 

- [SJ T 

AF at + AR at 

AF at + n AF bt 

AF bn - AR m 

J3 

< 


AF bt - AR at 
AF b 


( 11 ) 


( 12 ) 
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Premultiplying equation (9) by [S 1 ] T and substituting from equation (11) gives the reduced system 


AF m + AF t 


AF at + AF bt 


- |S,] T (K][S 2 | - [SHIRKS, 


Closed and sliding nodes (tangential movement) . — In this case equations (5) are constraining relations 
that need to be satisfied. After substitution from equations (5), equations (8) take the form 


AF an + AR an 


AF at ± n AR a 
AF bn - AR a 
AF bt =F MAR a 


A ^an + 


By following a similar procedure to that of case (2), equations (14) can be reduced to the form 


AF m + AF t 


AF at aAF m - [Q 3 ][K ][Q 2 ] “ = [Q 3 ][K HQ,] AU at 
AF bt ± /*AF an „ AU bt 
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where [Qj], [Q 2 ], and [Q 3 ] are transformation matrices defined as 
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( 16 ) 


Observe that in this case the stiffness matrix [Q 3 ][k][Qj is not symmetric. The dissymmetry is due 

to the Coulomb friction forces. The ± and signs are chosen to dissipate energy because the work due 
to friction forces is negative. 


The governing equations, equations (8), (13), and (15), can be solved for the incremental displace- 
ments {AU^}, {AU at }, {AU b }, and {AU bt }. Then by backsubstitution the incremental forces 
{AR^}, {AR at }, {AR b }, and {AR bt } can be obtained. However, at each incremental load the assumed 
contact conditions must be checked and adjusted, if necessary. Table I lists the contact region criteria 
and the transition boundaries. Iterations must be performed at each increment until convergence between 
assumed and calculated contact and frictional conditions is obtained. 


AUTOMATION OF METHOD 

The analysis method is based upon the assumption that the size of the load increments is sufficiently 
small that there is a linear load-displacement relationship and also that there is no more than one change 
of contact condition or phase in any load step. An algorithm to automatically limit the size of the load 
increment can be developed by studying the transitions between contact phases. The procedure is based 
upon an evaluation of the amount of load increment needed to reach a phase boundary. This “amount” 
is measured in terms of “scale factors.” 

The development proceeds as follows: First, we categorize the contact phase change as being either 
(1) from open to contact or (2) from contact to open. Next, within the contact region we categorize 
phase change as being either (3) from sticking to sliding or (4) from sliding to sticking. 

From open to contact . — Suppose a load increment {AF} causes an open pair of contact nodes, say a 
and b, to come into contact and even pass each other (analytically). Then, due to the assumed linearity, 
the proportion a^AF} of the load increment needed to exactly close the gap distance 6 can be deter- 
mined. Specifically, let the relative displacement of a and b in the direction n, normal to the surface, 
due to {AF} 1 be AU^ — AU bn . Then the scale factor needed to close the gap without penetra- 
tion is determined by the relation 

“i(AU m - AU bn ) = 6 (17a) 
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or 


“1 - «/(AU„ - AU bn ) 


(17b) 


From contact to open . — If the normal force R n between two contacting nodes a and b becomes a 
tensile force, the surfaces will separate. If AR n is the normal load increment during separation, in a 
typical iteration step, the proportion a 2 AR n of the load increment needed to exactly reach the transi- 
tion between contact and separation is given by the relation 

R„ + a 2 AR n = 0 or a 2 = _R n/ AR n ( 18 ) 


where a 2 is the scale factor. 

From sticking to sliding contact . — From Coulomb’s law of sliding with friction there will be a transi- 
tion between sticking and sliding at a contacting node pair when the tangential force R t exceeds /xR n , 
where is the friction coefficient and where, as before, R n is the normal force. Let AR t and AR n 
be increments in the tangential and normal contact forces due to a load increment {AF}. If the load 
increment is sufficiently small, AR t and AR n are proportional to the magnitude of {AF} and hence 
proportional to each other. Then the proportions c* 3 AR t and a 3 AR n of the tangential and normal 
force increments needed to reach the transition between sticking and sliding are given by the expression 

±(R t + a 3 AR t ) = fi( R n + a 3 AR n ) (19) 


where the sign is positive if R t AR n < R n AR t and negative otherwise. Hence, the scale factor a 3 is 


and 


a 3 


R t “ ^ AR n 

AR t - \i AR n 


for R t AR n < R n AR t 


( 20 ) 


a 3 


R t + P AR n , 

for 

AR t + /x AR n 


R t AR n ^ R n AR t 


( 21 ) 


Note that the static (stick condition) friction coefficient is generally slightly higher than the dynamic 
(sliding) friction coefficient. This difference has been neglected in this analysis. 

From sliding to stick contact . — It is difficult to establish and calculate a load scale factor for the 
transition from sliding to sticking contact. Therefore, for this case we simply assume all sliding nodes are 
sticking as a new load increment is applied. Then, we examine the sliding/sticking character of the con- 
tacting nodes: If they are sliding, no status change has occurred. However, if they are sticking (indi- 
cating a status change), the load increment must be reduced to establish the location and loading of the 
phase change. 

After the scale factors are determined for each case, for all contacting node pairs, the mini- 

mum of these (designated a) is selected as the global load scale factor. That is, the load increment 
{AF} is reduced to «{AF} and the process is repeated. 
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When friction is present and where there is a tentative phase change from sticking to sliding, the 
accuracy of the incremental procedure can be improved by monitoring the friction forces R t . Specifically, 
sliding will occur when R t approaches jiR n . Hence, an additional iteration criterion can be established 
by requiring R t to satisfy the relation 

(1 - e)/iR n < R t < fiR n (0 < e < 1) (22) 


where e is an accuracy parameter. 

Figure 3 shows a flow chart outlining the steps of the procedure. 


VALIDATION: CONTACT BETWEEN IDENTICAL INFINITE CYLINDERS 

A computer code was written to execute the foregoing algorithm. To validate the algorithm and the 
code, we first considered the contact between two semicylinders as depicted in figure 4. The semicylin- 
ders are geometrically identical, but they were given different elastic properties to investigate and illu- 
strate the effect of elasticity upon the contact mechanics. They each have radius R a = R^ = 1.0 inch 
and infinite length. Their elastic moduli and Poisson ratios are 

E a = 3.0xl0 7 psi (steel) 

E b = 1.8x10 psi (beryllium copper) 

j/ a = 0.292 ( 23 ) 

= 0.285 


The cylinders were aligned with parallel axes and pressed together by a uniform distributed load p of 
5x10 psi, as depicted in figure 4. Figure 5 shows the finite element mesh for the semicylinders. The 
mesh is very fine in the contact region. It has a total of 412 nodes and 364 linear quadrilateral plane 
strain elements. Three contact conditions were examined: Frictionless contact, friction contact without 
sliding, and friction contact with sliding. 

Frictionless contact . Figure 6 shows the resulting stress distribution for frictionless contact. The 
horizontal and vertical normal stresses and the maximum shear stresses are plotted along the vertical axis 
of symmetry. The stresses are measured in multiples of the maximum normal contact pressure p 0 of the 
Hertz theory. Comparisons with the classical Hertz solution (ref. 9) are also given. 

Figure 7 shows the distribution of normal stress in the contact area and also a comparison with the 
Hertz solution, where a is the semilength of contact of the Hertz solution. The contact area determined 
by the finite element solution is approximate in that the boundary between the contact and open regions 
lies between the closed and open node pairs. The boundary position can be estimated by interpretation. 

Table II presents a comparison of the numerical results and the Hertz results for the maximum con- 
tact pressure, the maximum shear stress, and the semilength of the contact. 
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Friction contact without sliding . — To determine the effects of friction, the stresses were calculated 
with friction coefficients /i of 0.0, 0.06, 0.12, and <». Figure 8 shows the results for the contact forces 
on the upper and lower cylinders. Figure 9 shows detailed representations for portions of the normal and 
tangential contact forces. Observe in figure 9(a) that the maximum normal contact forces increase with 
the friction coefficient. However, there is relatively little increase once the coefficient exceeds 0.12. 
Observe also that the normal contact forces decrease slightly in the edge region when the friction coef- 
ficient increases. 

Finally, table III shows the results for the semilength of the contact. As expected, the length 
decreases as the friction coefficient increases. 

Friction contact with sliding . — To determine the effect of friction when there is slipping or sliding 
between contacting nodes, let the upper cylinder of figure 4 rotate counterclockwise relative to the lower 
cylinder. Figures 10 and 11 show the resulting maximum compressive (principal) stresses and maximum 
shear stresses on the upper cylinder in the contact area for various friction coefficients. Analogous and 
similar results were obtained for the lower cylinder. 


CONTACT STRESSES BETWEEN SPUR GEAR TEETH 

Consider a pair of identical spur gear teeth in contact with each other at their pitch points. Let the 
teeth be modeled as in figure 12 with the tooth on the right being the driving tooth. Let the teeth have 
a fillet radius of 0.045 inch, a face width of 0.25 inch, and a pressure angle of 20°. Let the teeth be ele- 
ments of gears with a diametrical pitch of 8.0 and with 28 teeth per gear. Let the elastic modulus be 
30xl0 6 psi with a Poisson’s ratio of 0.3. Let the driving gear have a clockwise torque of 50 000 in. -lb. 

The stresses were calculated for the same friction cases as before. 

Frictionless contact . — Figure 13 shows stress contour results for the maximum compressive (principal) 
stresses and for the maximum shear stresses. The maximum contact pressure was found to be 8.11 xlO 5 
psi, which is approximately 10 percent higher than that predicted by the Hertz method (ref. 10): 

Friction contact without sliding . — Figure 14 shows the distribution of normal and tangential contact 
forces for the driving gear for four friction coefficients: 0, 0.05, 0.10, and Analogous results were 
obtained for the driven gear. 

Friction contact with sliding . — Figure 15 shows the maximum compressive (principal) stress and 
maximum shear stress distributions in the contact region for the driving gear tooth. Analogous results 
were obtained for the driven gear. 


DISCUSSION 

The results demonstrate the feasibility and the practicality of using the finite element method for gear 
stress calculations. Success in obtaining reliable results, however, is dependent upon the incrementation 
of the load to establish the extent of the contact area, which is unknown at the start of the problem. 

Because the contact region boundaries are unknown, they must be determined iteratively. These 
boundaries are more sensitive to the loading than are the stresses. This is seen in table II, where the dif- 
ference in the stresses between the Hertz and finite element analyses is much less than the difference in 
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the contact region. The sensitivity of the contact region arises from the near-parallel surfaces of the con- 
tacting bodies in the contact region. 

The inclusion of friction effects creates an additional complexity in the analysis. Indeed, friction pro- 
duces a dissymmetry in the stiffness matrix that significantly increases the computational effort. The 
computational effort was reduced by treating the friction forces as known external tangential forces that 
are determined from the previous iteration step. 

The results also show that friction produces higher stress in the contact region than when friction is 
absent. This means that precise surface geometry, polished surfaces, and lubrication can significantly 
reduce the stresses and thus increase surface life. 

The gear stress analysis is not comprehensive but only representative and illustrative of spur gear 
tooth stresses. The analysis is conducted for pitch point contact where there is no sliding due to gear 
kinematics. This means that an assumption of frictionless contact at the pitch point is very reasonable. 
Indeed, figures 14 and 15 show that friction has a relatively small effect on the analysis. 

Plane strain finite elements were used for the analysis presented in this paper. Plane strain is most 
suitable for very thick bodies such as the infinitely long cylinders used to compare the analysis with Hertz 
theory. For gear teeth, plane stress finite elements may be more appropriate, especially for narrow-face- 
width gears. Finally, the analysis assumes either static or sliding contact, whereas the contact between 
gear teeth in service is divided between sliding and rolling phases. Hence, more study is needed to obtain 
a comprehensive stress analysis. 


CONCLUDING REMARKS 

An incremental finite element procedure has been investigated for the analysis of contact stresses with 
application to meshing spur gears. The following conclusions have been reached: 

1. The feasibility and practicality of the procedure were established through several examples and by 
favorable comparison of results with those of the classical Hertz method. 

2. Because the contact area is not known a priori, it was necessary to employ an iterative procedure 
(herein called the load incremental procedure) to simultaneously determine the stresses and the contact 
area. 

3. The presence of friction between contacting surfaces increases the contact stress. This in turn 
means that precise surface geometry, polished surfaces, and effective lubrication can reduce contact 
stresses and thus increase surface life. 

4. The method is directly applicable to gear stress calculations. However, more analyses are needed to 
obtain a comprehensive understanding of contact mechanics throughout a mesh cycle. 
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TABLE I.— CONTACT CRITERIA AND 
TRANSITION BOUNDARIES 


Condition 

Criterion 

Before 

After 

Open 

Open 

(AU bn - AU.J - f n > 0 ! 

Close 

(AU bn - AU an ) - « n < 0 

Stick 

Stick 

Kn > o. |RJ < HR.„I 

Open 

R an<° 

Sliding 

Kn > l R .t! > M|R an | 

Sliding 

Stick 

Ran > 0- AR at (AU at - AU bt ) > 0 

Open 

R an< 0 

Sliding 

Ran > 0. ARa,(AU at ~ AU b J < 0 


TABLE II.— COMPARISON OF RESULTS OF FEM AND HERTZ 



FEM 

Herts 

method 

Error, 

percent 

Maximum contact 
pressure, p 0 , ksi 

19.33 

19.76 

-2.20 

Maximum shear stress, 

r max» 

5.536 

5.928 

-2.60 

Semilength of contact, 
a, in. 

.01732 

.01611 

7.51 


TABLE III.— FRICTION 
COEFFICIENT AND 
SEMILENGTH OF 
CONTACT 


Friction 

Semilength 

coefficient, 

of contact, 


a, 


in. 

0 

0.01732 

.06 

.01730 

.12 

.01729 

00 

.01727 
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Figure 3. — Flow chart for load incrementation iteration. 
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Rgure 4. — Contact between two semicylinders. 



Figure 6. — Subsurface stresses along axis of symmetry. (Stress 
values are normalized by p 0 , the maximum normal contact 
pressure.) 


lalized stress, P(xVp 0 



Figure 7. — Distribution of norma! stress in contact surface. 
(Stress values are normalized by p 0 , the maximum normal 
contact pressure.) 


Nodal contact forces, lb 


Coefficient 
of friction, 
P 


0 

O .06 

□ .12 




Figure 8.— Distribution of contact forces. 
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of friction, 



Tangential distance from Initial contact point along 
contact area, in. 

Figure 10. — Maximum compressive (principal) stresses on 
contact surface of upper cylinder (sliding case). 



(b) Tangential contact forces (upper cylinder) — an 
enlargement of area II of figure 8(a). 


Figure 9.— Detailed representation of forces. 
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Figure 1 1 . — Maximum shear stresses on upper cylinder in con- 
tact area (sliding case). 
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Figure 1 2.—Spur gear teeth In contact at pitch point. 



(a) Maximum compressive (principal) stress. 



(b) Maximum shear stress. 

Figure 1 3. — Stress contours for Motionless contact. 


19 


Contact forces, (b 









REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, 
gathenng and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this 
collection of information, including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington. VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188). Washington, DC 20503. 


1. AGENCY USE ONLY {Leave blank ) 


4. TITLE AND SUBTITLE 


2. REPORT DATE 

March 1992 


Contact Stresses in Meshing Spur Gear Teeth: Use of an 
Incremental Finite Element Procedure 


6. AUTHOR(S) 

Chih-Ming Hsieh, Ronald L. Huston, and Fred B. Oswald 


3. REPORT TYPE AND DATES COVERED 

Technical Memorandum 


5. FUNDING NUMBERS 


WU-505-63-36 
1L16221 1A47A 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 


E-6776 


10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 


NASA TM- 105388 
AVSCOM TR-90-C-029 


7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

NASA Lewis Research Center 
Cleveland, Ohio 44135-3191 
and 

Propulsion Directorate 

U.S. Army Aviation Systems Command 

Cleveland, Ohio 44135-3191 


9. SPONSORING/MONITORING AGENCY NAMES(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 

Washington, D.C. 20546-0001 

and 

U.S. Army Aviation Systems Command 
St. Louis, Mo. 63120-1798 


11. SUPPLEMENTARY NOTES 

Chih>Ming Hsieh and Ronald L. Huston, Department of Mechanical, Industrial, and Nuclear Engineering, University 
of Cincinnati, Cincinnati, Ohio 45221-0072 (work funded under NASA Grant NSG3-188). Fred B. Oswald, NASA 
Lewis Research Center. Responsible person, Fred B. Oswald, (216) 433-3957. 


12a. DISTRIBUTION/AVAILABIL1TY STATEMENT 

Unclassified - Unlimited 
Subject Category 37 


13. ABSTRACT (Maximum 200 words) 

Contact stresses in meshing spur gear teeth are examined. The analysis is based upon an incremental finite element 
procedure that simultaneously determines the stresses in and the contact region between the meshing teeth. The 
teeth themselves are modeled by two-dimensional plane strain elements. Friction effects are included, with the 
friction forces assumed to obey Coulomb’s law. The analysis also assumes that the displacements are small and that 
the tooth materials are linearly elastic. The analysis procedure is validated by comparing its results with those for 
the classical two contacting semicylinders obtained from the Hertz method. Agreement is excellent. 



14. SUBJECT TERMS 


Contact stress; Finite element; Spur gears; Friction 


15. NUMBER OF PAGES 

22 


16. PRICE CODE 


17. SECURITY CLASSIFICATION 18. SECURITY CLASSIFICATION 19. SECURITY CLASSIFICATION 
OF REPORT OF THIS PAGE OF ABSTRACT 

Unclassified Unclassified Unclassified 


20. LIMITATION OF ABSTRACT 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. Z39-18 
298-102 





















